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04 . Hadron spectroscopy has always played an important role in the revealing mechanisms underlying the dynamic of 
^""^ ■ strong interactions. 

\ The heavy hadron containing a single heavy quark is particularly interesting. The light degrees of freedom (quarks, 
^ ^ . antiquarks and gluons) circle around the nearby static heavy quark. Such a system behaves as the QCD analog of 
O ' familar hydrogen bound by the electromagnetic interaction. 

2^ In Refs. [l|,l4| relativistic generalization of the three-body Faddeev equations was obtained in the form of dispersion 

relations in the pair energy of two interacting quarks. The mass spectrum of S'-wave baryons including u, d, s 
^f) , quarks was calculated by a method based on isolating the leading singularities in the amplitude. We searched for the 
approximate solution of integral three-quark equations by taking into account two-particle and triangle singularities, 
all the weaker ones being neglected. If we considered such an approximation, which corresponds to taking into account 
two-body and triangle singularities, and defined all the smooth functions of the subenergy variables (as compared 
1 ' with the singular part of the amplitude) in the middle point of the physical region of Dalitz-plot, then the problem 
Q-l was reduced to the one of solving a system of simple algebraic equations. 

' In Ref. 3 the relativistic six-quark equations are found in the framework of coupled-channel formalism. The 
dynamical mixing between the subamplitudes of hexaquark are considered. The six-quark amplitudes of dibaryons 
are calculated. The poles of these amplitudes determine the masses of dibaryons. We calculated the contribution of 
six-quark subamplitudes to the hexaquark amplitudes. 

In the present paper the six-quark equations for the nonstrange baryonia with the open charm are found. The 
nonstrange baryonia BBc are constructed without the mixing of the quarks and antiquarks. The six-quark amplitudes 
. of baryonia are constructed. The relativistic six-quark equations are obtained in the form of the dispersion relations 

■ over the two-body subenergy. The approximate solutions of these equations using the method based on the extraction 
. of leading singularities of the amplitude are obtained. The paper is devoted to the calculation results for the baryonia 

■ mass spectrum (Table |T|. In conclusion, the status of the considered model is discussed. 

] The model in question has only two parameters of previous paper gluon coupling constant go = 0.314 and cutoff 
parameter Aq = 11. We used the cutoff Age = 5.18 which is determined by AI ~ 4100 MeT^ (the threshold is equal to 
4130 Mey). 

The quark masses of the model are = 495 MeV and rric = 1655 MeV. The estimation of theoretical error on 
the S'-wave hexaquarks masses is 1 AleV . This results was obtained by the choice of model parameters. 
, We consider 9 baryonia with the content qqQqqq and the spin-parities = 0~, 1~, 2~ . The isospins are equal to 

■ 1 ) f J I (Table HI . We have predicted the masses of baryonia containing c quark using the coupled-channel formalism. 
We believe that the prediction for the S'-wave charmed baryonia is based on the relativistic kinematics and dynamics 
which allow us to take into account the relativistic corrections. 

The quark pairs Qq use there not so many, therefore the baryonium masses cannot increase enough with the 
decreasing of the cutoff Age. 

The degeneration of baryonium masses with the different spin-parities = 0^, 1^ was obtained. We cannot also 
calculate the bound states of baryonia with = 3^. 

The baryonium state [uucddd) for the spin-parities = 0^, 1^, 2^ is calculated with the nine subamplitudes: 
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seven ai (similar to a\ ) and two a2 ^ , a" ^ 
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The baryonium (uucudd) consists of 16 subamplitudes with the spin-parities J ==0 ,1 ; 12 ai and 4 a2. 
a\ ^ , Q.\ ° , a2 ^ J Q^2 ° ■ For the case of the spin-parity ~ 2^ the subamplitude ° is absent. The 
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TABLE I: qqQqqq, q = u,d, Q = c. 
Quark masses rUq = 495 MeV, nic = 



Parameters of model: cutoff A = 11.0, Age = 5.18, gluon coupling constant g = 0.314. 
1655 MeV. 
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EcA, AEc, KN, NK, KA, AE*, AcA, AAc 
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states with spin-parities J =0 ,1 ,2 (udcuuu) are constructed with 13 subamphtudes: 10 ai and 3 q;2: a2 
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. The baryonium udcuud for the spin-parities J = ,1 takes into account 23 subamphtudes: 17 

^Quc-^uu^ ^Qdc-^uu ^j^^ ^^^^ _ 2_ ^j^^ subamphtudes a2"'''^"'') 



ai and 6 ^2; 0"" -0"°°"'" -O^^o 
Q!2 ° , ° are absent. 

The system of equations of the baryonium SgA {uuc ddd) for the spin-parity = 1^ (as the example) was 
constructed: 
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af " = A + 2 af° 7i(l""0"^) + 6 af '7i(l"''l'"^") + 6 af'7i(P"0'"^') , 

= A ar"/i(0"n"") + af'^/i(0"'=0"'=) 3ar''/i(0""r'^") 3 af" /i(0""0"'^") 

+ 3 '/i (0"^r^") + 3 af'/i (0"^0^'*") , 

= A + 2 "li (I'^^'^^l'^"'^") + 4 af'/i (I'^^'^^l"'^") + 4 af" Ji (I'^^'^^O"'^") -h 2 "ji (l'^^'^^^'^") 

= A + af" 7i (l"**"!"") + af°/i (1"'*'0"'=) + 2 ""/i (l"'*"l<*"'*') + 3 al""" hirh""^) 

+ 3af '/i(l"V'^") + ar'/i(r'^"r'^") + af ' hir'^O^'^) + 2 a^"^^'' hil'^h'^^l'^'^) 

+ 2ar'''-f2(l"''"0"'=l'*"'), 

= A + af" Ji (0"-^^!"") + af 7i (0"'^"o"") + 2 a}"" Ji (0"'^"l'^"'^") + 3 of ' Ji (0"'^"l"'*") 

+ 3 af'h (O^dQud) al"" li{0''h'=^) + /i(0"'^"o"'^") + 2 ai""!"" J2(0'"^"l'"'l'^''^') 

= A 2a}'"/i(l'^'^"l'^"'^") + 2ar'hil''h^'') 2 af" /i(l'=V'^") -h 2af'/i(l'='^"r'^") 
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2afhil^'^0"^) +4af^^''l2{l"^0"n^^) , (6) 
A + 2ar^i(0'='^l'*'^) + 2af '/i(0"'*"l'"^") + 2 af" /i(0'='*>'*") + 2 al"" hiO'^h''^) 
2af''li{0''^0''^) +Aa°"''^'"'l2{0''^0''n^^) , (7) 
A + 2af''/4(r"l'^'^0"'=) + 2al'"'h{r''l'^'^l'^'^) + 4 ar'/3(l""l'^"'^"l'"^") 

4 0°"" /3(1""1'*%"'^") + 4af'i'"/6(l™'l'^"'*"0'"=l'^"'^") , (8) 
A + af "/4(0"'=l'^"'^"l™) + /4(0""1'^"'^>") + 2 a}'"'/4(l'^^"'l'^"''") 
2af'' /3(0'"=l'^"'^"l"'^") + 2 0?"" /3(0'"=1'*"'*>'*") + 2al'"'l3{0'"'l^h''^) 

2 af" 73(0"^'*^'='*") + 2 al'^"^'"'le{0"n^h''"V'^^) + 2 a°2"'"^'"' hiO'^'V^^O"''!^^) . (9) 
We used the functions Ji, I2, I3, I4, Iq similar to the paper Q. 

The poles of the reduced amplitudes a; correspond to the bound states and determine the masses of the charmed 
baryonia. 

In Table H] the calculated masses of nonstrange baryonia with the open charm are shown. 

We predict the mass of lowest charmed baryonium with the isospin I = \ and the spin-parity = 1^ [M = 
3293 Me V^). 

The known way with which to calculate the low-energy properties of hadronic systems rigorously is Lattice QCD 
(LQCD) 0, @|- In LQCD calculations, the quark and gluon fields are defined on a discretized space-time of finite 
volume of the lattice volume, such deviation can be systematically removed by reducing the lattice spacing, increasing 
the lattice volume and extrapolating to the continuum and infinite volume limits using the known dependences 
determined with effective field theory (EFT) 0-01 ■ 

We try to consider the tasks which are similar to the Lattice calculations. 
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